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Abstract 

A micromechanical formulation is presented for the prediction of the overall thermo- 
inelastic behavior of multiphase composites which consist of short fibers. The analysis is 
an extension of the generalized method of cells that was previously derived for inelastic 
composites with continuous fibers, and the reliability of which was critically examined 
in several situations. The resulting three-dimensional formulation is extremely general, 
wherein the analysis of thermo-inelastic composites with continuous fibers as well as 
particulate and porous inelastic materials are merely special cases. 


1 Introduction 

Modeling of the constitutive response of advanced composite materials continues to be an 
important part of these strategic materials’ development. Analytical models that predict the 
effective behavior of composites are used not only by mechanicians in structural analysis of 
large-scale composites but also by materials scientists in developing new material systems. 
For an analytical model to fulfill these two distinct functions, it must be based on the 
micromechanics approach which allows one to generate the average response of a composite 
material from the properties of the individual constituents and their geometric arrangement. 
Only then can such a model be used by a material scientist to investigate the effect of different 
deformation mechanisms on the overall response of the composite in order to identify the 
appropriate constituents for a given application. At the same time, a micromechanical model, 
if it is to be used in a large-scale structural analysis program, must be computationally 
efficient in addition to being able to generate accurate displacement and stress fields at 
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both the macro- and micro-level. Knowledge of accurate field quantities at the micro-level 
is important in predicting failure of a structural element. 

A number of models presently exist that can fulfill one or the other of the aforementioned 
tasks. On the other hand, there are very few working models that are both computationally 
efficient and sufficiently accurate at the micro- as well as the macro-level to satisfy the distinct 
needs of the material science and structural mechanics communities. One such model is the 
method of cells (Aboudi (1991)) and its generalization (Paley and Aboudi (1992), Aboudi 
(1993)) which is believed to have the potential of fulfilling both aforementioned tasks. 

In the original formulation of the method of cells, a continuously-reinforced, unidirectional 
fibrous composite is modeled as a rectangular, doubly-periodic array of fibers embedded in 
a matrix phase. The periodic character of the assemblage allows one to identify a repeating 
unit cell that can be used as a building block to construct the entire composite. The proper- 
ties of this repeating cell are thus representative of the properties of the entire assemblage. 
The unit cell consists of a single fiber subcell surrounded by three matrix subcells. Hence 
the name method of cells. The rectangular geometry of the repeating unit cell allows 
one to obtain an approximate solution for the stresses and strains in the individual subcells 
given some macroscopically homogeneous state of strain or stress applied to the composite. 
The approximate solution to the thus posed boundary- value problem is, in turn, used to 
determined macroscopic (average) or effective properties of the composite. The macroscopic 
behavior is displayed in terms of: effective elastic moduli; effective coefficients of thermal ex- 
pansion; effective thermal conductivities; and effective stress-strain response in the inelastic 
region. 

In the generalized method of cells for continuous fibrous composites, the repeating unit 
cell can consist of an arbitrary number of phases. Hence the generalized method of cells is 
capable of modeling a multiphase composite. This is of importance since even a two-phase 
composite with a metallic constituent behaves after the development of plasticity effects as 
a multiphase composite. A user’s guide for a computer program based on the generalized 
method of cells has been recently presented by Aboudi and Pindera (1992). 

In the original formulation of the method of cells for short-fiber composites, the unit cell 
consists of a single region occupied by the inclusion, surrounded by seven distinct matrix 
regions. 

In the present paper the generalized method of cells is developed and formulated for 
the modeling of multiphase thermo-inelastic short-fiber composites. The presented three- 
dimensional micromechanical analysis consists essentially of four steps. It starts by iden- 
tifying a repeating volume element of the periodic multiphase composite, followed by the 
definition of macroscopic average stresses and strains from the microscopic ones. In the third 
step the continuity of tractions and displacements are imposed at the interfaces between the 
constituents. These establish, in conjunction with micro-equilibrium, the relationship be- 
tween microscopic total, thermal and plastic strains and macroscopic strains via the relevant 
concentration tensors. In the final step the overall macroscopic constitutive equations of 
the composite are determined. These four steps form the basis of micro- to- macromechanics 
analyses which describe the behavior of heterogeneous media (Suquet, 1985). The resulting 
three-dimensional micromechanical analysis establishes the overall elasto-plastic behavior of 
the multiphase thermo-inelastic composites. This is expressed as a constitutive relation be- 
tween the average stress, strain, thermal and plastic strains, in conjunction with the effective 
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elastic stiffness tensor. 

The micromechanically established constitutive equations are consistent in the sense that 
when they are employed to predict the effective coefficients of thermal expansion of the 
multiphase composite, they provide the same estimate obtained from Levin’s result (1967). In 
the latter, the overall coefficient of thermal expansion is expressed in terms of the properties 
of the constituents in conjunction with the mechanical concentration factors. Similarly, 
the derived constitutive equations are also consistent in the sense that the resulting overall 
inelastic strain of the multiphase composite coincides with the expression obtained by Dvorak 
(1992). Here the overall inelastic strain is expressed in terms of the inelastic strains of the 
phases in conjunction with the mechanical concentration factors. Both Levin’s and Dvorak’s 
results follow from the relationship between the overall and local transformation fields which 
can be established by invoking the reciprocal theorem (Dvorak and Benveniste (1992)). 


2 Model Description 


Consider a composite material with a periodic structure whose repeating volume element 
consists of Na x Np X jV 7 rectangular parallelepiped subcells. The volume of each one of 
the subcells is d a hpi yi) where a, ft and 7 are running indices: a = 1 , . . . , N a ; ft = 1, • • - Np 
; 7 = l,...,iV 7 , in the X\ , rr 2 and - directions, respectively. The total volume of the 

repeating volume element is dhl where d = d a , h = Ylp=i hp , t = Fig. 

1 an example is shown for a repeating volume element with N a = 3, Np = 4 and jV 7 = 2. 

Each of the subcells can be filled, in general, by an elastic-viscoplastic temperature - 
dependent material. Consequently, the overall behavior of the composite should simulate a 
multiphase thermo-viscoplastic composite. By appropriately selecting the subcells dimen- 
sions or, alternatively, by appropriately filling the subcells, various types of composites can be 
considered. Thus, unidirectional long-fiber composites, short-fiber composites, porous ma- 
terials and laminated materials can be modeled as special cases. Consequently, the present 
micromechanical multiphase composite model is very general and various situations can be 
obtained as special cases. 

By an approximate micromechanical analysis of the detailed interactions of the subcells 
of the repeating volume element, overall constitutive relations which govern the effective 
behavior of the multiphase viscoelastic composite, can be established. This analysis relies 
on the requirements that static equilibrium of the materials in the subcells are ensured, 
and that continuity of the displacements and tractions between neighboring subcells within 
the repeating volume element, as well as between neighboring repeating volume elements is 
satisfied on an average basis. 

Let us introduce local coordinates x^\ £3^) whose origin is located at the center 

of the subcell (a/37). These local coordinates are shown in Fig. 2 for subcell {ocft 7) and the 
neighboring one (<*7/?) in the Xi - direction, where a is defined to have the form 


f a + 1 a < N a 

\ 1 a = N Q 


( 1 ) 


This definition ensures that for a < N a the neighboring subcell in the x\ - direction is 
the one labelled by (a + l,/?,7) within the repeating volume element, whereas for a = N a 
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the neighboring subcell is within the next repeating volume element whose first subcell is 
(l/?7). Similarly, /3 and 7 are defined by 


0 


-(?*' 


0<N fi 

P = Np 

(2) 

II A 

-i -j 

(3) 


Since the average behavior of the composite is sought, it is sufficient to consider a first 
order theory in which the displacements in the subcell are expanded linearly in terms 

of the distances from the center of the subcell, i.e., in terms of x^\ and x^* Thus, the 
following first order expansion in the subcell (a/?7) is considered 


u\ aplf) = + x[ a) <f>^ ] + X { 2 0) x\ a ^ ] + * = 1,2,3 


(4) 


where w\ a ^ are the displacement components at the center of the subcell, and , X% a ^ 

and are microvariables that characterize the linear dependence of the displacement 

u\ aPy) on the local coordinates e( l n « (4) and the sequel, repeated Greek 

letters do not imply summation. Note that due to the linearity of eqn. (4) , static equilibrium 
of the material within the subcell (a/3 7) is ensured. 

The components of the small strain tensor are given by 


4f 7) = + dju^) 

Ld 


ij = 1,2,3 


(5) 


where d\ = 8/8x [ a ^ , 82 — d/dx ^ and 83 = 8/8x 3^ . 

The average strains in the subcell, are given according to eqns. (4)-(5) by 


~(a0-i) _ i(o,0-y) ~(a0-y) _ (a0y) ~(a0- y) _ j(a0-y) 

£ 11 — 9\ 1 e 22 — A2 1 6 33 — Y3 

2e<f 7) = xi alh) + 4 a/h) 

2 eSf 7) = V’i° /5 ' y) + 4 0/?7) 


2e<f 7) = X [ alh) + <A ath) 


(6) 


The average strains in the composite are expressed as 


N a Nfi TV, 


C *i ~ S 2C dahpty 

unc o=l 0=1 7=1 


(7) 


The constitutive equation of the elastic- viscoplastic temperature-dependent material that 
occupies subcell (a/3 7) is expressed as a relationship between the averages of stress 

total strain ej“^ 7 \ inelastic strain ejj a ^\ and the thermal strain = aj“^A T ( where 

j s coefficient of thermal expansion tensor, and AT is the temperature deviation 
T — Tr from a reference temperature Tr) , i.e. 


-(a/0-y) _ ^,(a0'r)/-(a0~t) -I(a 0y) -T(a0'i)\ 

a ij ~ U ijkt \ e ki ~ e kt ~ e kt ) 


( 8 ) 
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where is the elastic stiffness tensor of the material. 

The average stress in the composite is determined from 


(T,j — 


dht 


N a N e N 1 

a—l p=\ 7=1 


> -(or^-y) 

7 °ij 


( 9 ) 


It will be shown in the following that by employing the displacement and traction conti- 
nuity conditions at the interfaces between the subcells of the repeating volume element, and 
at the interfaces between neighboring repeating volume elements, it is possible to eliminate 
the micro variables and obtain, via a smoothing operation (homogenization), a set of con- 
tinuum equations that model the overall behavior of the multiphase short-fiber composite. 
This is achieved by establishing relationships which connect the microstrains at the subcells 
to the total overall macrostrains in the composite via the appropriate concentration tensors. 


3 Interfacial Continuity of Displacements 

At any instant, the displacement components must be continuous at the various interfaces 
of the repeating unit cell, and at the interfaces between neighboring repeating cells. This 
implies that for a = 1 , . . . , N a ; 0 = 1 , . . . , Np ; 7 = U • • • » , the following relations hold 


u 


(apy) 


X 


(<*) _ ^ 


U 


(a/?7) 


-(<*) _ 


X\ ' = 


( 10 ) 


u 


(apy) 


*(« - bn 

X 9 — 


U 


(apy) 


*(/*) _ 
X 9 


-/i- 


( 11 ) 


,(« 07 ) 


-U) _ ix U i 

X Q O 


(aPy) 


X 3 ” 2 


( 12 ) 


Notice that with the definition of a, jd and 7 , eqns.(l)- (3), continuity of the displacement 
at the interfaces between neighboring repeating cells is ensured. 

Continuity conditions (10) - (12) are imposed at the interfaces in an average sense. For 


example, eqn. ( 10 ) is applied 

in the 

form 


M 2 

rl-t/2 1" 


- 

J — hpl*! * 

l~ty/ 2 [ 


~(«) _ rfn 

x i - 2 J 

y/i/ 3/2 

ft- r /2 

(cr/?7) 


J-hp/2 

J-t-y/ 2 

u] 

*(“) _ 
x \ — 2 


dxi^dx^ = 


dx-f^dx^ 


Using eqn. (4) in (10) one obtains 


w\ apy) + \d Q ^ ] = - \d~<t>f M) 

It £ 


(13) 


(14) 
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( 15 ) 


Two similar equations would result from the use of eqn. (4) in (11) and (12): 


1 W*’ = 

2 2 


(16) 


All field variables in eqn. (14) are evaluated at the center line Xj°^ of the subcell (a/3 7), 

and at the center line ip of the subcell (a/3 7), see Fig. 2. In the final form of the displace- 
ment continuity relations it is required to evaluate all terms at a common location. To this 
end, let x^ denote the location of the interface between subcell (a/3 7) and the following one 
(a/?7). It follows that 


r <0 - -d 

1 2 a 


(17) 


In addition, the location of the center line Xj of the subcell ( S/? 7) with respect to x^ 
is given by 


xf 0 = x</> + |<i s 


(18) 


By expanding the field variables in eqn. (14) in Taylor’s series about the common interface 
x</> and dropping second order terms, we obtain 


w\ a0l) - -d a (^-w { ^ l] - 4 aPl) ) = + ^(iT-w^ - #^>) 


d 


(19) 


2 ~'dx x * ^ ' ‘ 2 ay dx t 

where all field variables in (19) are evaluated at the common interface x^ . The discrete 
structure of the composite is eliminated by a smoothing operation according to which eqn. 
(19) is assumed to be valid at any point of the homogenized medium, which effectively 
represents the multiphase composite, for all a, /? and 7. For a composite subjected to 
homogeneous boundary conditions, the behavior of all repeating cells is identical and a 
uniform field exists at the equivalent homogeneous medium. 

Let F/ a ) be defined as: 


where 


F} a) = w\ a0y) + f! a) - wf M) + ff a) 


fi a) = - *r* 7) ) 


2~“ v dx 1 ' 

Consequently, equations (19) can be rewritten as 


F (a) = 0 


a — 1 , . . . , N a 


(20) 


( 21 ) 


( 22 ) 


A similar derivation based on eqns. (15) -(16) would respectively provide the following 
two relations 
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G\ p) = 0 , fi = l ,...,Np (23) 

and 

W,W=0 , 7=1,...,JV, (24) 

where the following definitions have been employed 

Gf> = w\ a/h) + g\ p) - w\ ah) + gf (25) 

H \ 7) = w\ af3y) + - w\ a ^ + fcp ( 26 ) 

s!'” = -\ h ^ w '?' h) - x!" M ) < 27 ) 

W = - # W ) (28) 

Z C/^3 

Like eqn. (22), eqns. (23)-(24) are simultaneously valid at any point of the homogenized 
medium, which effectively represents the multiphase composite subjected to homogeneous 
boundary conditions, for all a, /? and 7. 

From eqns. (22) - (24) we obtain, respectively, that 

= 0 , 53Gjfl= 0 , S^ ) = 0 ( 29 ) 

cr=l /?=1 7=1 

The three summations provide, respectively, that 


~ r 

£/!“> = <> . £*!« = o . £^’ = 0 

cv— 1 /?=1 7=1 


Since 


• a^’ = 0 (31) 

for all a , /? and 7 , it follows by differentiation of eqns. (22), (23) and (24) with respect to 
Xi , X 2 and £3, respectively, that 


d ( 
~Z~ W i 

OX 1 


^ (or)J 


(a/?7) ^ (o/?7) 


^ (a/?7) ^ (of/?7) 

OX 2 5^2 

d ...(ag-r) _ ^ 

d* 3 *' ' 
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Equations (31) are satisfied by assuming that common displacement functions, to,-, exist 
such that 


_ w . (33) 

for all a, /? and 7. This assumption is consistent with the fact that the entire repeating 
volume element of the periodic composite is mapped by a smoothing operation into a single 
point within the equivalent continuum medium. The displacement components at this point 
are u? t . 

Using eqns. (33), we readily obtain from (30) the following set of continuum relations 


N 3 


E = * 


0 =\ 


N-, 


dx 5 


E wl *™ = 

7=1 


-Wi 


dx 3 


P=l,...,Np 

7=1,...,^ 

(34) 

II II 

H- 1 I-* 

(35) 

a = 1, . . . , N a 

P = 

(36) 

terms of the microvariables <f>\ al3 ' r \ X; 0 ^, 


It is possible to derive an equivalent system of equations which are expressed in 
terms of the average subcell strains 

The average composite strains, e,j, were defined in eqn. (7). Let us first show that 


_ 1 ( dwi t dwj x 

= o( q_ + Q_ ) 


(37) 


2 y dxj dxi 

For i = j = 1, let us multiply eqn. (34) by and perform a summation over /? from 
1 to N f 3 , and over 7 from 1 to N y . This provides 


No, N 0 ^7 O,,, 

E E E = dkt? w ' 

a=l ( 3 = 1 7=1 


dxi 


(38) 


Comparing eqn. (38) with (7) and using the first equality in (6), gives In = 

For i = 1, j = 2 let us multiply eqn. (34) with i = 2 by hp( y and perform a summation 
over (3 from 1 to Np, and over 7 from 1 to N y . Similarly, let us multiply eqn. (35) with 
i = 1 by d a £ y and perform summations over all a and 7. By adding the resulting relations 
we obtain 


E + x[° M ) = + 

Of*/?,' 7 1 



(39) 


Again, using the last equality in (6) in conjunction with equation (7) and comparing this 
to the above expression gives €12 = (|r*- + 1^-)/2. In a similar manar the other four relations 
in (37) can be established. 
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It is now possible to express the continuum eqns. (34)-(36) in terms of and e,j. 

Setting * = 1,2 and 3 in (34)-(36), respectively, we obtain 


E = A,, , 

a — 1 

II II 

H-i t— ‘ 

(40) 

E b l 22^ = hl 2 2 , 
0=1 

P 

II II 
J”" ‘ 

(41) 

Ny 

2 > ^7^33 — ^-^33 ? 

7=1 

a = 1, . . . , N a 
ft = 1,..., Np 

(42) 

The addition of eqn. (34) (with i = 2) multiplied by hp and summed 
(with i = l) multiplied by d a and summed over all or, yields 

over all /?, to (35) 

Na N P 

E E + xi“ w 

a=l p—1 



i.e., 



N a Up 

E E d a hp^ = dht\i 

a=\ P-l 

II 

(43) 

which are the desired continuum equations expressed in terms of subcell and average strains. 
Similar operations give 

E E 

0=1 7=1 

, a = 1, . . . , N a 

(44) 

Na Ny 

E E = die 13 

a=l 7=1 

£ 

II 

<50. 

(45) 


Equations (40) - (45), form a set of N a Np + NpN^ + N a N 7 + N a + Np + relations 
which replace eqns. (34)-(36). These relations can be written in a matrix form as follows 


A g e, = Je (46) 

where the 6-order average strain-rate vector is defined by 

c = (cil 5 ^22 1 £33 , 2e23 ) 2ci3 , 2 E 12 ) (47) 

and the 6N a NpNy order subcell strain vector, € s , is defined by 

€ a = (6< 1U \...,6<^ W ’>) , (48) 

where the 6 components of the vector are arranged as in eqn. (47). The matrix Ag is 

N a (Np+N 1 +l)-\-Np(N 1 + \)+N^ by 6N a NpN y , while J is a N a (Np+N y +l)+Np(N^+l)+N y 
by 6 matrix. It should be noted that the matrix Ag involves the geometrical dimensions of 
the repeating cell only. 
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4 Interfacial Continuity of Tractions 

The tractions must be continuous at the interfaces between the subcells of the repeating 
cells, and at the interfaces between neighboring repeating cells. These conditions, imposed 
in an average sense, give 

-(aP'l) _ 
a \i ~ a U 


~{a0i) _ -(or/J 7 ) 


>2 i 


= 02 , 


(49) 


- (qt/3-y) 


' 3, 


= a 


(apy) 


3, 


where i — 1,2,3 ; a = 1 ; 0 = 1 and 7 — l,...,iV 7 . It can be easily 

verified that these equations involve some repetitions. It can be shown that the system of 
independent interfacial conditions is 


=.(*07) — 
a n — a u 

«=!,-• 

.,N a - 

1 



P = 





7 = 1 ,.. 

-,N^ 


(50) 

- (or/37) _ 
a 22 — a 22 

« = !,•• 

.,N a 




/?=!,-• 

■,N 0 - 

1 



7 = 1 ,.. 

■ ,N . r 


(51) 

— (or/? 7 ) _ 
a 33 — °33 

«=!»•■ 

..,iVa 




/3 = 1 ,., 





7 = 1 ,.. 


1 

(52) 

- (a/?7) _ -(a0 7 ) 
a 23 “ °23 

a = 1 , . 

..,iv a 




/?= 1 ,. 


1 



7 = 1 ,.. 

,.,N y 


(53) 

-(otpi) _ -(a Pi) 
a 32 — °32 

a = 1 ,. 

..,N a 




P = N 0 

7 = 1 ,- 

■ ■ ,N~, — 

1 

(54) 
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(or/hf) __ 

13 “ °13 


a = 1 , . . . , N a — 1 

0 = 1, 

7 = 1,...,AT 7 


(or/57) _ 

31 — °31 


a = N a 

0 = 1 , 

7 = 1,..., AT 7 - 1 


( 55 ) 


(56) 


r 12 


= a 


(a/37) 


12 


a = 

0 = 1,...,^ 

7 = 


(57) 


-(a/?7) _ -(or/57) 


'21 


21 


a = N a 

0=1 ,...,W/5-l 

7 = 1, . . . , ATy 


(58) 


As in the case of the displacement continuity conditions, the continuum model requires 
that for a composite subjected to homogeneous boundary conditions, conditions (50) - (58) 
are satisfied at all points of the homogeneous equivalent medium. By using the constitutive 
laws (8), it is possible to represent eqns. (50)-(58) in matrix form, that is, 




(59) 


where the 6N a NpN y - (N a Np + N a N^ + NpN^) - (N a + N 0 + AT 7 ) by 6N a NpN^ matrix A M 
involves the elastic properties of material in the subcell, and 


= (e / ( 111 ),. 

g I(JV«JV,JV 7 )) 

(60) 

= (€ T < m >,. 

t f{N a N p N y )^ 

(61) 


5 Overall Inelastic Thermomechanical Constitutive Law 

The combination of equation (59) and (46) leads to 

Ae, - £>(ef + ej) = Ke (62) 

where 
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(63) 


0 

J 

Solving eqn. (62) for the subcell strains e s yields 

e, = Ae + D(e[ + ej) 

where 


A = 


Am 

A g 


D = 


Am 

0 


K 


A = A X K , D = A 1 D 

Let the mechanical concentration matrix A be partitioned into N a NpN~y 6-order square 
submatrices in the form 


Similarly, let 


A = 


A (m) 


D = 


D< m > 1 


D (N a N^) 


where £)^ a ^ are square matrices of 6 N a Np N y order. 
It follows from eqn. (63) that 


(64) 


(65) 


+ D {a<h \e\ + ej) (66) 

Equation (66) expresses the average strain in the subcell in terms of the uniform overall 
strain e (the applied macrostrain) and the subcell inelastic and thermal strains, via the con- 
centration matrices A^ a01 ^ and £f( a0 ' , \ Notice that the same matrices are operating 

on both the inelastic and thermal strains of the subcells. 

It was shown by Dvorak and Benveniste (1992) and Dvorak (1992) that for any rep- 
resentative volume element under a uniform overall strain e and temperature change AT, 
which contains a piecewise uniform distribution of thermal and inelastic fields (eigenstrains) 
associated with the applied loading e , the averages of the local strain can be expressed in 
terms of the mechanical and eigenstrain concentration tensors. These tensors depend on the 
local elastic moduli, and on the shape and volume fraction of the phases and are therefore 
constant. Equation (66) is consistent with this representation of the average subcell strain, 
where the constant mechanical and eigenstrain concentration tensors are given by A^ a0 ^ 
and D^ a0ri \ respectively. 

Substitution of equation (66) into (8) yields 

0 .(<*P'r) = 7 )(ef -f ej) — (c / (°'^'d -(- (67) 

Consequently, in conjunction with eqn. (9), the following effective elastoplastic thermo- 
mechancial law of the composite can be established 
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( 68 ) 


& = B*{e - e 1 - e T ) 

where the effective elastic stiffness tensor, B *, of the composite is given by 

1 Na N, , Ny 

B- = Tu 7 E E E 

ant o=l /3=1 7=1 

and the composite inelastic strain tensor is defined as 


(69) 


e' = 


-B 


,_1 N a N $ N- 


dhl 


E E E d a h^C^\D^e{ - e 1 ^) 


ar=l 0=1 7=1 

and the average thermal strain tensor as 


(70) 


D—l JV„ N-, 

g g 2 d a h^C^\D^e T s - e T ^>)) 


£ = 


dht ^ ^ 

uuc a=l 0=1 7=1 


(71) 


The effective coefficient of thermal expansion vector, or* of the composite, is given by 


a" = e r /A T (72) 

It can be numerically verified that the effective thermal stress vector r* = B*ot * is given 
in accordance with Levin (1967) result by 

= 317 E E E <h„t,A (73) 

“h" a=i/?=i7=i 

where A T ^ a ^ is the transpose of A^ a ^\ and is the thermal stress vector of the 

material filling subcell (a/? 7 ) . This implies that the generalized method of cells model 
is consistent. This consistency follows from the fact that when this model is employed to 
predict the effective coefficients of thermal expansion, it provides the same estimate obtained 
from Levin’s result (73). 

Like eqn. (73), it is possible (e.g. Dvorak (1992)) to establish the following estimate for 
the overall inelastic strain e 1 of the composite 


1 


7V q Nfi N~f 


i'-asEEE-t.W 

unc 0=1/3=17=1 


,7(007)^/(007) 


(74) 


where Q T (°'^'d is the transpose of the mechanical concentration tensors, that connect 
the stress in the subcell to the total overall stress & of the composite subjected to 

homogeneous boundary conditions, i.e. 


It can be easily verified that 


_ Q(<*0 7)^. 


(75) 


Q(or07) __ 5* -1 


(76) 
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As in the thermal case, it can be numerically verified that the generalized cells model is 
consistent in the sense that the overall inelastic strain determined directly from eqn. (70), 
coincides with the values obtained by employing eqn. (74). 

The micromechanically established overall thermo-inelastic constitutive law (68) is valid 
for any type of thermomechanical loading (i.e., any combination of normal, shear and thermal 
loadings). A significant advantage of this constitutive law stems from the fact that it does 
not rely on any symmetry conditions that may exist under certain types of applied loadings. 
Thus in the implementation of this law, the question as to whether such symmetry conditions 
exist or not is irrelevant. 


6 Conclusion 

A micromechanical theory is offered according to which an elastic-thermo-inelastic consti- 
tutive law which governs the overall behavior of short-fiber composites is established. This 
theory relies on the assumption that the composite possesses a periodic structure so that 
it is sufficient to analyze a repeating volume element. This repeating volume element (or 
cell) can include an arbitrary number of subcells, so that multiphased composites contain- 
ing short fibers and arranged in various configurations (e.g. square, hexagonal, or square 
diagonal packing) can be considered. 

The established constitutive law is based on the principles of continuum mechanics where 
microequilibrium and continuity of displacements and tractions at the various interfaces are 
ensured in an averaged sense. The key ingredient in the construction of this macro constitu- 
tive law, is the establishment of the appropriate concentration tensors, A and at 

the microlevel. Given these concentration tensors, the average stress and strain micro field 
variables in the subcell can then be expressed in terms of the uniform overall strain and the 
subcell inelastic and thermal strains, see equations (66) and (67). The macromechanical law 
is then established by applying a homogenization (or averaging) procedure, see equations (7) 
and (9), to the microscopic thermo-inelastic constitutive equations. 

The established analytic constitutive law can be readily applied to investigate the behav- 
ior of various types of composites given knowledge of the behavior of the individual phases. 
Numerous advantages can be stated regarding the current macro/micro constitutive laws as 
compared to the other numerical micromechanical approaches in the literature, e.g. the finite 
element unit cell approach. One advantage is the fact that any type of simple or combined 
loading (multiaxial state of stress) can be applied irrespective of whether symmetry exist or 
not, as well as without resorting to different boundary condition application strategies as in 
the case of the finite element unit cell procedure. Another, advantage concerns the availabil- 
ity of an analytical expression representing the macro elastic- thermo -inelastic constitutive 
law, thus ensuring an economy of memory when implementing this formulation into a struc- 
tural finite element analysis code. Furthermore, this formulation has been shown to predict 
accurate macro behavior given only a few number of subcells, within the repeating cell (see 
Paley and Aboudi (1991), and Arnold et al.(1993)). Whereas, if one employs the finite ele- 
ment unit cell procedure, a significant number of finite elements are required within a given 
repeating cell to obtain the same level of accuracy as with the present formulation. As a re- 
sult it is possible to utilize the presented constitutive law to efficiently analyze metal matrix 
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composite structures subjected to complex thermomechanical loadings. This is particularly 
important when analyzing realistic structural components, since different loading conditions 
exist throughout the structure, thus necessitating the application of the macromechanical 
equations repeatedly at these locations. 
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Fig. 1 . (a) A repeating volume element that consists of N a = 3, Np = 4 and N-y = 2 
subcells. 
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